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ABSTRACT 

The effects of ambipolar diffusion on the linear stability of weakly ionised accretion discs 
are examined. Earlier work on this topic has focused on axial magnetic fields and perturbation 
wavenumbers. We consider here more general field and wavenumber geometries, and find that 
qualitatively new results are obtained. Provided a radial wavenumber and azimuthal field are 
present along with their axial counterparts, ambipolar diffusion will always be destabilising, 
with unstable local modes appearing at well-defined wavenumber bands. The wavenumber 
corresponding to the maximum growth rate need not, in general, lie along the vertical axis. 
Growth rates become small relative to the local angular velocity when the ion-neutral col- 
lision time exceeds the orbital time. In common with Hall electromotive forces, ambipolar 
diffusion destabilises both positive and negative angular velocity gradients. In at least some 
cases, therefore, uniformly rotating molecular cloud cores may reflect the marginally stable 
state of the ambipolar magnetorotational instability. 
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1 INTRODUCTION 

The behaviour of the magnetorotational instability (MRI) in 
weakly ionised gas is a problem of considerable importance, bear- 
ing directly on our understanding of the dynamics of protostel- 
lar discs. The low ionisation regime is characterized not only 
by increased Ohmic dissipative losses (Jin 1996), but by Hall 
(Balbus & Terauem 2001; Wardle 1999) and ambipolar diffusion 
>Blaes & Balbus ..19 94) effects as well. The latter two effects are 
of particular interest, since they can be destabilising under certain 
circumstances. 

The Hall and ambipolar diffusion regimes are generally im- 
portant in distinct regions of parameter space ( Balbus & Terauem 
I2OOII) : the former is dominant (along with Ohmic losses) at the high 
densities typical of inner regions of protostellar discs, whereas the 
latter dominates in the outer regions of such discs (where they are 
typically observed), and under interstellar conditions. In this paper, 
except for a brief excursion in the Appendix, we confine ourselves 
to the ambipolar regime, and study the general axisymmetric be- 
haviour of the MRI in these environments. 

The linear and nonlinear effects of ambipola r diffusion on 



ideal and Hall MHD, it turns out that it is not generally true in the 
ambipolar limit. Here, the fastest growing wavenumbers may have 
both vertical and radial components, and it is necessary to perform 
a full axisymmetric analysis to identify these modes. 

The plan of the paper is as follows: In §2, we discuss the basic 
formulation of the problem, including a presentation of tempera- 
ture and number density regimes associated with non-ideal MHD 
effects. Section 3 is a linear analysis of differentially rotating discs 
for axial wavenumber and field configurations. This serves to unite 
this paper with previous work on ambipolar diffusion. Section 4 
extends this analysis for arbitrary field geometries and perturbation 
wavenumbers, and contains the key results of the paper. In §5, we 
present a brief summary. 



2 PRELIMINARIES 

The fundamental fluid equations are mass conservation, 



the MRI have been studied in both the two-fluid ^lae^^^albu^ 
19941 : lHawleY& Stone I Il998h and single-flui d jMac Low et~ 



IhtawleviK Stonelliwai) ana smgte-llui a IMac Low et al. 
1 9951 IWardle 1 1 999t ISalmeron & Wardle |2003^ regimes. Our pur- 



dp 
dt 



+ V-(p«) = 0, 



pose in revisiting this problem is to highlight a novel result of prac- 
tical interest. Previous studies were confined to vertical wavenum- 
bers only, based on the assumption that this direction corresponds 
to the fastest growing local modes. Although this is true for both 
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and the equation of motion. 



(1) 

B. 

(2) 



Our notation is standard: p is the mass density, v the velocity, P 
the gas pressure, B the magnetic field, and <E> the gravitational po- 
tential. The density, velocity, and pressure all refer to the dominant 
neutral species. 
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Figure 1. Regions of ambipolar (A), Hall {H), and Ohmic (O) dominance . 
Numbers are taken from eqs. (24) and (25) fronj| Balbus & Terauemli200lh . 
assuming v\/cg = 0.1. Note that there is no explicit dependence on radius, 
central mass, or ionisation fraction. The region of ambipolar dominance of 
interest is the left-most wedge. 



dt 



The induction equation is (see e.g., Balbus & Terquem 2001) 
47rryJ JxB {JxB)xB 



Vx 



vxB 



(3) 



Here 77 is the resistivity, c the speed of light, J the current den- 
sity, Tie the electron density, e the magnitude of the electron 
charge, 7 the electron-ion drag coefficient, and Pi the ion mass 
density. Numerical value s of rj and 7 are (.Draine et al...l983^ 
iBalbus & Terquem l200lh : 



■n = 234 



T 



1/2 2 -1 
' cm s , 



: 2.75 X 10"cm^s"^g"\ 



where n is the number density of the neutrals and T is the temper- 
ature. Finally, for future reference, we define the Alfven velocity in 
the usual way: 



(4) 



The terms on the right side of equation ^3) correspond respec- 
tively to Faraday induction, Ohmic resistivity, Hall electromotive 
forces (HEMFs), and ambipolar diff usion. The relative ratio o f the 
ambipolar to Hall terms is given by iBalbus & Terquem 1200 ll) : 
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(5) 



where Ca is the isothermal sound speed. Assuming that the final 
two factors are each about 0.1, we see that a neutral density below 
about 10'' cm~^ brings us into the ambipolar diffusion regime (see 
Fig. 1), and we shall assume that this restriction is satisfied. We 
retain the resistivity in §3, however, to illustrate the simple relation 
of ambipolar diffusion to Ohmic dissipation for the case of axial 
field geometry, and in the Appendix, where both Hall and Ohmic 
terms are included in the derivation of a very general dispersion 
relation. 



3 AXIAL FIELDS AND WAVENUMBERS 
3.1 Stability 

We consider first the local stability of a differentially rotating disc 
threaded by a weak vertical field, B — Bez,m the presence of am- 
bipolar diffusion. We use standard cylindrical coordinates {R, 0, 2) 
with the origin at the disc centre, and consider plane wave pertur- 
bations that depend only on z. Linearized quantities (indicated by S 
notation) are proportional to exp{at — ikz), where k is the vertical 
wave number. In the Boussinesq limit, this corresponds to fluid dis- 
placements in the plane of the disc, so vertical structure is unimpor- 
tant. Under these circumstances, pressure, density, vertical velocity, 
and vertical magnetic field perturbations all vani sh. The two-fluid 
versio n of this problem has been considered bv iBlaes & Balbusi 
il994l) : we shall compare the results of our single-fluid treatment 
with those of these authors. 

Stability is most easily assessed by working in the limit a 
0, assuming that the stable-unstable transition is bridged by a cr = 
solution. The linearized equations of motion are 



ikB 
ikB 



SBr = 0, 



0, 



where k is the epicyclic frequency defined by 

2 ,o2 , dn^ 

The linearized induction equations are 

k^ri'SBR - ikBSvr = 0, 



k -q'SB^ 



where 



dlnR 



SBn — ikBSv^ = 0, 



V =v-^ 

is the Ohmic resistivity modified by ambipolar diffusion. 
Solving the above system of equations, we find 

2,2/2, 2 f dQ,^ , , 2 2 A „ 



(6) 
(7) 

(8) 

(9) 
(10) 

(11) 
(12) 



This is the desired marginal stability condition. It is easily shown 
that the left hand side should be negative for instability. The insta- 
bility criterion is then, 



,22^ dn 



din 7? 



1 + 



2/2 1 
K 1] 



(13) 



a more restrictive condition than the classical MRI. 

It is often the case that Ohmic resistivity is small com- 
pared with ambipolar d iffusion, which was the limit assumed by 
lBlaes& Balbus H1994 . If we take 77 ^ 0, 77' wi/7Pi, the in- 
stability criterion becomes 



,22^ dii 
k vai < 



dlnR \ + y^Pi 
where we have introduced the ion Alfven velocity 



VAi 



VA 



1/2 



(14) 



(15) 



The analogous two-fluid instability criterion obtained by 
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iBlaes & Baibus1<1994 was 



A;^«A, < — - — =- I — I fBB94 instability criterion) 

dlnR \ K.^ + 72pi J J J 

(16) 

Clearly, the inequality <14t is consistent with the inequality <16t 
provided that 

IP > 



(This ignores the effect of ambipolar diffusion in the unperturbed 
background disc.) Since explicitly enters the dispersion relation 
and we are assuming a simple exponential time behaviour for the 
perturbations, we set Ba = in the analysis for self-consistency. 
The inclusion of a Br component should not qualitatively affect 
the stability behaviour, though it may affect growth rates. We shall 
also ignore Ohmic losses, as ambipolar diffusion generally domi- 
nates over this term (cf. Fig. 1). 



a quantitative implementation of the assumption of negligible ion 
inertia. 



4.2 Dispersion Relation 

The linearized mass conservation equation in the Boussinesq limit 
is: 



3.2 Axial Wavenumber Dispersion Relation 

We next obtain the full dispersion relation for finite a. The lin- 
earized radial and azimuthal equations of motion are now 



ikB 
ikB 

4:TVp 



SBr = 0, 
SB^ = 0. 



The same components of the linearized induction equation are 



(ct + k'^r]')SBR - ikBSvr = 0, 



(a + k^r]')SB, 



4> 



SBr - ikBSz 



0. 



dlnR 

The resulting dispersion relation is 

a* + 2k'^rj'a^ + Sacr^ + 2k^rj'{K^ + k'^vl)a + Bo = 0, 
where the constants B2 and Bo are given by 

dO.'^ " 



B2 



2 I ni 2 2 , I 4 /2 , 
K + 2k Vj^ + k Tj + 



dlnR 



Bo = k^v% 



dln_R 



,22 

+ k va 



, 2,4/2 

+ K k r; . 



(17) 
(18) 

(19) 
(20) 

(21) 

(22) 
(23) 



Ambipolar diffusion acts on axial wavenumber disturbances only 
as though it were a field-dependent additive resistivity, ri' . Adding 
other field components while retaining k — kez, 01 adding other 
wavenumber components while retaining B = Bez, does not 
change the dispersion relation. Qualitative changes are introduced 
only if nonaxial wavenumbers as well as nonaxial field components 
are considered. 



4 GENERAL ANALYSIS 
4.1 Preliminary Assumptions 

We now consider the axisymmetric behaviour of the instability with 
more general field geometries and wavenumbers. We ignore buoy- 
ancy, so our analysis holds either either for a barotropic disc, or 
locally at the midplane. The perturbation wavevector is 



k = kRBR + kzBz, 



(24) 



kRSvR + kzSvz = 0. 



(26) 



Using this to eliminate Svz from the equations of motion leads to 
the radial and azimuthal forms 



k k-B 

<t5vr — -j^2Q.5v4, + "I^^SBr = 0, 



47rp 



K k-B 

<^Sv^ + -tttSvr, — i- — SBa, = 0. 
20, Anp 



(27) 



(28) 



By similarly eliminating SBz via a vanishing divergence condi- 
tion, only the R- and ^-components of the induction equation are 
needed. The linearized induction equations are 



^Tvy pip 
kRB^ik-B 



iiv^Pip 



SBr - i{k-B)5vR 

SB4, = 0, 



SB^ - i{k-B)St 



dO, 
dlnR 



+ 



kRB^{k-B) 



SBr = 0. 



The resulting dispersion relation is 

cr'' + Cacr^ + C2cr^ + Cia + Co = 0, 



with 



Ca = 



C2 = 



Ci = 



fc^B^ + (k-B) 
4-Kjpip 



I^Z 2 



2{k-VA) + 

kRB^{k-B) dVL 



kB(k-B) 



4:11-1 Pi p 



4:TT-ypip dlnR 
(k-VA) + -p-zt 



C3 



Co = (k-VA) 



{k-VA) + ^2 

fc| .2 



kl dO? 1 



, i^z 2 



{k-VA) 



A:2 dlnR 
kRB^{k-B) dO. 



(29) 



(30) 



(31) 



(32) 



(33) 
(34) 



kB{k-B) 



4-K"ipiP 



ATYjpip dlnR 



(35) 



and disturbances have space-time dependence exp(ife-r + at). In 
the presence of shear, a toroidal magnetic field B^ would grow 
linearly with time if a radial field component Br were also present: 

dO, 



B^it) = B4,{0)+tBR 



dlnR' 



(25) 



4.3 Stability 

We assume that the transition from stability to instability proceeds 
through a — 0, a point that can be confirmed numerically. In the 
limit (T ^ 0, 



4 MM Kunz & S.A. Balbus 



a = -Co/Ci. (36) 
With Ci > 0, the instabihty requhement is then 

Co < 0, 

From equation <35> . it is apparent that Co consists of the usual MRI 
terms plus explicit ambipolar diffusion terms. The first of these is 
simply an effective resistivity, but the second introduces a novel 
form of coupling: in formal Cartesian index notation k) it is 

kiBjdvj /dxi. 

It is the dominant ambipolar term in Co when 7Pi is large, and can 
always be chosen to be destabilising. Since this term is nonvanish- 
ing only in the presence of a nonaxial wavenumber component, it 
is easy to show that the fastest growing wavenumber need not, in 
general, lie along the 2-axis. 

This has important consequences. Let us view the marginal 
stability transition Co = as an equation for 



X 



{k-VA.) 
f22 



as a function of kn/kz, for a given set of disc parameters. Transi- 
tion to instability is present only if there are solutions to this equa- 
tion with X real and positive. With 



one finds that the Co 



X 



fe| dlnn 



~2 2 ir\2 

It — Hi /\l , 

requirement is satisfied when 

2 _ — — — 

-ypi kz Bz 



> n 



kn B^ dlnfi 
kz Bz d\nR 



The denominator vanishes at a wavenumber ratio 



kz 



n B^ dlnO, 
'ypi Bz dluR 



1 



' Bz 'ypi 



(37) 



(38) 



It is straightforward to show that there will always be solutions 
with large positive values of X near this selected value of kn/kz- 
These large wavenumber modes are well localised in a WKB sense. 
In other words, ambipolar diffusion always destabilises differential 
rotation, whether the profile is increasing inward or outward. De- 
spite the presence of a resistive term in the dispersion relation, there 
are unstable modes at large nonaxial wavenumbers. Moreover, be- 
cause fluid motions do not simply bend field lines locally when 
ambipolar effects dominate (loss of flux-freezing), some of the sta- 
bilising effects of large wavenumber magnetic tension seen in the 
standard MRI are lost. 

The modification of the MRI in the presence of ambipolar dif- 
fusion does not change the fundamental cause of the instability, 
which is based on rotational dynamics: angular momentum is re- 
moved from fluid elements with less angular momentum and given 
to fluid elements with more angular momentum, an intrinsically 
destabilising process. What is complicated by the presence of am- 
bipolar diffusion is the fluid element "tethering," now no longer a 
simple, spring-like magnetic tension based on flux-freezing. The 
radial component of the disturbed magnetic field, for example, de- 
pends both upon azimuthal as well as radial motions, and field lines 
do not follow fluid elements. One consequence is the destabili- 
sation of outwardly increasing rotation profiles. In astrophysical 
settings, the primary interest is not outwardly increasing angular 
velocity profiles of course, but the additional destabilisation that 




Figure 2. Instability growth rates in the (k-VA.)^ /fl^ (ordinate), \k/kz\ 
(abscissa) plane for a Keplerian disc with B^/Bz = 5 and Q/fpi = 
0.1, 0.5, and 1 (from top to bottom). Only regions of instability are shown, 
with the height proportional to the growth rate. The maximum growth rates, 
from top to bottom respectively, are 0.410, 0.170, and O.llO. Note that 
the classical MRI maximum growth rate 0.75f2 occurs at (k-VA.)'^ /Q^ = 
15/16 andfc/fez = 1. 



attends the usual outwardly decreasing profiles. Ambipolar diffu- 
sion ^oesjiotjjiowever, lead to growth rates in excess of 0.75Q, 
lBal bus&Hawlevlll992h . Indeed, as shown in the following sec- 
tion, growth rates become a small fraction of fl when the ambipolar 
parameter i^/fpi significantly exceeds unity. 

4.4 Growth Rates 

It is of interest to calculate explicitly some representative growth 
rates in the two-dimensional phase space described by the 
wavenumber parameters (k-VA)/^ and k/kz- In Figure 2, we 
present three-dimensional plots of the growth rate for representa- 
tive cases with fl/"/pi ^ 1. While the general tendency of ambipo- 
lar diffusion is to lower the maximum growth rate below 0.75^, 
there is also a significant extension of the domain of unstable non- 
axial wavenumbers. 

As Q/jpi approaches and exceeds unity, the most unstable 
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Figure 3. As in Fig. 2, but for VL/^pi = 2, 5, and 10 (from top to bottom). 
The maximum growth rates, top to bottom respectively, are 0.08f2, O.OIQ, 
and 9 X IQ-^n. 

modes are associated with a near cancellation of the resistive and 
nonresistive (destabilising) ambipolar terms in the Co term of the 
dispersion relation. The modes appear at progressively larger val- 
ues of kn/kz- The maximum growth rate decreases sharply from 
~ O.IQ. at SI ~ 7p, to ~ O.OOm at ~ 107Pi (see Fig. 3). 
Once the rotation frequency drops below the ion-neutral collision 
frequency, the growth time becomes very long compared with an 
orbital time. The instability remains viable, however, if its growth 
rate remains shorter than other evolutionary timescales. 



5 SUMMARY 

In this paper we have examined the axisymmetric magnetorota- 
tional instability in the presence of ambipolar diffusion. Our results 
should find their primary applicability in molecular discs on in- 
terstellar or galactic scales, as well as the low density outermost 
regions of protostellar discs. We have allowed radial and axial 
wavenumber geometries, and azimuthal and axial magnetic field 
components. The inclusion of a radial magnetic field component 
would induce a linear time dependence in the azimuthal field com- 



ponent. Since appears explicitly in the dispersion relation, we 
considered only the case of a vanishing background radial field 
component. 

The combination of nonaxial wavenumbers and magnetic 
fields greatly enhances the response of a differentially rotating gas 
to ambipolar diffusion, compared with the axial geometries that 
are usually analysed. In particular, the wavevector corresponding to 
the maximum growth rate is nonaxial when ambipolar diffusion is 
present. Short wavelength unstable wavenumber modes are always 
present, even for increasing outward angular velocity profiles. In a 
non-accreting system, solid body rotation is an energy extremum. 
Differential rotation causes no local linear instabilities in a purely 
hydrodynamic gas, and even in the presence of a magnetic field, 
ideal MHD only destabilises outwardly decreasing angular velocity 
profiles. It is therefore noteworthy that non-ideal MHD processes 
open paths to destabilise any differentially rotating configuration. 
A likely consequence of the action of the instability in nonaccret- 
ing systems is solid body rotatio n, which may explain the velo city 
profiles of molecular cloud cores iSarranco & Goodman Il998h . In 
a disc in which accretion is a possibility, the nonlinear outcome of 
the instability will not be solid body rotation, but pronounced dif- 
ferential rotation and turbulence. 

The parameter that measures the relative importance of am- 
bipolar diffusion is Q./'ypi, the ratio of the angular rotation veloc- 
ity to ion-neutral collision rate. Although there are unstable large 
wavenumber modes at all values of this ratio, the growth rate of 
the instability drops rapidly below Q, when this ratio approaches 
or exceeds unity. The viability of the instability then depends upon 
whether the growth time remains small compared to the expected 
lifetime of the system in question. When a background Bn is 
present, there is the additional uncertainty of whether the dynami- 
cal consequences of a steady, linear-in-time growth of B^ will out- 
pace a more leisurely growing exponential instability. 

A subsequent paper will address the stability of a system in 
which Ohmic losses, HEMFs, and ambipolar diffusion are simul- 
taneously present, the dispersion relation of which is given in the 
Appendix. 
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APPENDIX 

We present here the generalised dispersion relation including all 
non-ideal MHD effects: Ohmic diffusion. Hall effect, and ambipo- 
lar diffusion for the case Be. = 0, as in the text. Although the 
range of parameter space where all three terms are comparable is 
very small, it is nevertheless useful to have a general dispersion re- 
lation which may be conveniently specialised to a particular gas. 
Additionally, the final form of the complete dispersion relation has 
a mathematical form which can be written as a natural general- 
isation of the (comparativ ely) simple Hall dispersion relation of 
iBalbus & TerauemH200lh . This too is worth noting. 

Let us first introduce some preliminary notation. We start with 
the generalised resistivity parameters, ^'i and '3' 2: 
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(39) 



and their respective arithmetic and geometric means, T and ^f^: 

*i*2. (40) 



Next, there are two parameters with dimensions of magnetic 
field divided by frequency, /u, and B~ /v: 



V V Tie y tJcM B^pi 



(41) 



where tJc/j is the cyclotron frequency of a particle of mass ju, the 
mean mass per neutral particle, 

_ 

and cr is the unit vector in the radial direction. 

With these definitions, the generalised non-ideal MHD induc- 
tion equations may be written: 

[o + fc'^*!) &Br + tl^it^-J-SB^ - i{k-B)5vR = 0, (42) 



kRB{k-B+ 



ATTpiy 

dn kzB{k-B+) 



+ 



i-KpU 



dlnR Ampv 
UBz - i(k-B)Sv4, = 



SBf 



(43) 



(a + fc''*!) SBz - ) gg _ i(^k-B)5vz = 0. (44) 

The dynamical equations <27> and J28> remain unchanged. 

The non-ideal MHD dispersion relation may be written in a 
form similar to that of the axisymmetric dispers ion relation of a 
disc in the Hall regime jBalbus & Teraueml200il) : 
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[k-VA) 



«-Z 2 



a + ICo = 0, 
(45) 



with 

y- 4 2 , x2 , ,4,T,2 , kzB{k-B- 

/C2 = -p-K +2{k-VA) + K * + 



kzQ.B{k-B+) , fc| dO.^ 



1-KpV 



+ 



fc2 dlni? 



(46) 



{k-VAT + 



2 , kznB{k-B+) , fc| "I 



2TTpV 



+ 



fc2 dlni? 

(47) 



4^2 27rp!/ 

A detailed analysis of this equation is deferred to a later paper. 
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